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ABSTRACT 

We investigate the suppression of the baryon density fluctuations compared to the 
dark matter in the linear regime. Previous calculations predict that the suppression 
occurs up to a characteristic mass scale of ~ 10 6 M Q , which suggests that pressure has 
a central role in determining the properties of the first luminous objects at early times. 
We show that the expected characteristic mass scale is in fact substantially lower (by 
a factor of ~ 3-10, depending on redshift), and thus the effect of baryonic pressure 
on the formation of galaxies up to reionization is only moderate. This result is due to 
the influence on perturbation growth of the high pressure that prevailed in the period 
from cosmic recombination to z ~ 200, when the gas began to cool adiabatically and 
the pressure then dropped. At z ~ 10 the suppression of the baryon fluctuations is still 
sensitive to the history of pressure in this high-redshift era. We calculate the fraction 
of the cosmic gas that is in minihalos and find that it is substantially higher than 
would be expected with the previously-estimated characteristic mass. Expanding our 
investigation to the non-linear regime, we calculate in detail the spherical collapse of 
high-redshift objects in a ACDM universe. We include the gravitational contributions 
of the baryons and radiation and the memory of their kinematic coupling before re- 
combination. We use our results to predict a more accurate halo mass function as a 
function of redshift. 
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1 INTRODUCTION 

The detection of the cosmi c microwave backgro und (CMB) 
temperature anisotropics l|Bennett et al.l [l996) confirmed 
the notion that the present-day galaxies and large-scale 
structure (LSS) evolved from the primordial inhomo- 
geneities in the density distribution at very early times. 
After cosmic recombination, the gas decoupled from its me- 
chanical drag on the CMB, and the baryons subsequently be- 
gan to fall into the pre-existing gravitational potential wells 
of the dark matter. Regions that were denser than average 
collapsed and formed bound halos. First the smallest, least 
massive objects collapsed, and later, larger objects formed 
through a mixture of mergers and accretion. The formation 
and properties particularly of early galaxies at high redshift 
are being actively studied i n anticipation of many expected 
obser vational probes (e.g.. iBarkana fc Loebll200ll ; iBarkanal 
l2006h . 

A well-known solution for the collapse of a halo that 
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consists of dark matter on ly in an Einste i n de Sitter (EdS) 
universe was presented by iGunn fc Gottl (|1972| ). This solu- 
tion considers a spherical region initially with a small uni- 
form overdensity compared to the background universe. As 
the universe expands, the overdensity expands slower than 
the background until it reaches a maximum radius, turns 
around, and collapses. The critical overdensity, in the cor- 
responding linearly-extrapolated calculation, marks the col- 
lapse time of a dark matter halo in this case as 8 C — 1.686, 
a value that does not depend on the halo mass or collapse 
redshift. The mathematical solution gives a singularity as 
the final state, but physically we know that even a small 
initial asymmetry will make the object stabilize with a fi- 
nite size after reaching a virial equilibrium between motion 
and gravity. 

Extensive work has been done on spherical col- 
lapse models, especially models that include a cosmo- 
logical constant or a dark energy background (e.g., 
Lahav ct al. 1991; Deshingkar ct al 1l200ll; lLokas fc Hoffman 
200 it iHorellou fc Bergel 120051 : iMaor fc Lahavl 120051 : IWanj 
2005); in particular, the cosmological constant A changes the 
above value of overdensity (5 C ) by about 0.6%. In addition, 
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many numerical simulations of the formation of primordial 
objects at z ~ 20-30 have be en performed. However, the ear - 
liest stars formed at z ~ 65 ijNaoz. Noter fc Barkana 2006), 
and even for a halo that collapsed at z ~ 30, 5 must have 
started significantly non-linear (~ 9%) even for a simulation 
that begins as early as z ~ 600. 

With the ACDM cosmological parameters 
i|Spergel et all l2006h . the contribution of the photons 
to the expansion of the universe cannot be neglected 
when considering the form ation of the first objects 
|Naoz. Noter fc Barkanal I2006T I. Moreover, the baryons 
have a non-negligible contribution compared to the dark 
matter, and their different evolution must be included in 
the collapse process. 

When considering the formation and properties of 
the first luminous objects, we must investigate the rela- 
tion between the bar yon and the dark matter fluctuations. 
iGnedin fc Hull || 19981 ) defined a fiducial "filtering mass" that 
describes the highest mass at which the baryonic pressure 
still manage s to suppress the linear baryonic fluctuations 
significantly. IGnedin! l|2000h extended the usefulness of the 
filtering mass to the fully non-linear regime by showing that 
it is also related to another characteristic mass scale - the 
largest halo mass for which the gas content is substantially 
suppressed compared to the cosmic fra ction. As we show be- 
low, if we follow previous calculat ions (|Gnedin fc Huilfl998l ; 
iGnedinllioOol : IGnedin etaH 12003). we find a characteristic 
mass at high redshift of ~ 1O 6 M0, approximately constant 
at z > 60 and decreasing only slowly with time afterwards. 
This is somewhat larger than the mass scale of the first ob- 
jects and suggests a potent effect on the formation of the 
first objects. 

Here we present an improved calculation of the charac- 
teristic mass that is mainly based on the improved calcu- 
lation of the baryon density and temperatur e fluctuations 
that we presented in iNaoz fc Barkanal (|2005T l . We first re- 
view the basic equations of linear perturbation growth (Sec- 
tion [271]) • We then divide the power spectrum into several 
different ranges of scales that are associated with large-scale 
structure (Section l2.2l) . the filtering scale (Section 12. 3[) . and 
small scales (Section 12. 4|l . Note that we define the filtering 
mass with a different normalization than in previous works, 
as explained in Section 12.31 For completeness we compare 
our calculation to the older, inaccurate approximation of 
a spatially-uniform sound speed along with other approxi- 
mations (Sections 12.51 and 12. 6[) . We use our results for the 
filtering mass to estimate the gas fraction in minihalos (Sec- 
tion !2.7p . In Section[3]we calculate in detail the critical over- 
density for collapse of halos that form at very high redshifts, 
following the evolution of perturbations outside the horizon 
(Section l3.1|) and inside it (Section l3.2p . We also predict the 
halo abundance at different redshifts (Section 13. 3|) . Finally, 
we summarize and discuss our results in Section 4. 

Our calculations are made in a ACDM universe, in- 
cluding dark matter, baryons, radiation, and a cosmological 
constant. We assume cosmological parameters matching the 
three ye ar WMAP data tog ether with weak lensing obser- 
vations l|Spergel et al.ll2006n . i.e., fi m = 0.299, Ov = 0.74, 
Qi = 0.0478, h = 0.687, n = 0.953 and a 8 = 0.826. We also 
consider the effect of current uncertainties in the values of 
cosmological parameters on some of our results, by compar- 
ing to the results with a different cosmological parameter 



set specified bv lViel et all (|2006l ): 0. m = 0.253, fl A = 0.747, 
n b = 0.0425, h = 0.723, n = 0.957 and a s = 0.785. These 
parameters represent typical 1 -a errors, in terms o f the pa- 
rameter uncertainties given bv lSpergel et al.l l|2006l) . 



2 LINEAR GROWTH OF PERTURBATIONS 
2.1 The Basic Equations 

INaoz fc Barkanal |2005h showed that the baryonic sound 
speed varies spatially, so that the baryon temperature and 
density fluctuations must be tracked separately. Thus, the 
evolution of the linear density fluctuations of the dark mat- 
ter (5<jm) and the baryons (5b) is described by two coupled 
second-order differential equations: 

Sdm + 2-ff<5dm = — Hq — ^p- (fbSb + /dm^dm) , 

5 h +2H5 h = ^-Hq—^- (fb8h + /dm<5dm) n~^~ ($b + S 

2 a A a z fi 

where Q m is the present matter density as a fraction of the 
critical density, k is the comoving wavenumber, a is the scale 
factor, fi is the mean molecular weight, Ho marks the present 
value of the Hubble constant H, and T and St are the mean 
baryon temperature and its dimensionless fluctuation, re- 
spectively. These equations can be derived by linearizing the 
continuity, Euler, and Poisson equations. The baryon equa- 
tion includes a pressure term whose form comes from the 
equation of state of an ideal gas. Th e linear evolution of the 
temperature fluctuatio ns is given by (jBarkana fc Loebll2005l ; 
INaoz fc Barkanall2005l l 

where x e (t) is the electron fraction out of the total number 
density of gas particles at time t, 5 1 is the photon density 
fluctuation, t y = 8.55 x 10 _13 yr _1 , and T 7 and 8t 1 are the 
mean photon temperature and its dimensionless fluctuation, 
respectively. Equation ([2]) results from the first law of ther- 
modynamics, where in the post-recombination era before 
the formation of galaxies, the only external heating arises 
from Compton scattering of the remaining free electrons 
with CMB photons. The first term of equation comes 
from the adiabatic cooling or heating of the gas, while the 
second term is the result of the Compton interaction. Note 
that prior analyses (e.g., [Peebles 1980; iMa fc Bertschingerl 
1995) had assumed a spatially uniform speed of sound for 
the gas, but this assumption is inaccurate. 

It is often useful to express the evolution in terms of two 
linear combinations of the dark matter and baryon fluctua- 
tions. Defining 5tot = fbSb + /dm<5dm (in terms of the cosmic 
baryon and dark matter mass fractions fb and /dm) and 
A = 6b — Stat (following the derivation in iBarkana fc Loebl 
2005), and using equations |T} we can write two differential 
equations that describe the evolution of S to t and A: 

5tot + 2HStat = t:Hq — ^-<5 t ot — h^— - — (<5tot + A + <5t) , 
2 a d a 2 (i 

A + 2HA = -/ dm ^^_(5 tot + A + «5 T ) . (3) 
Note that the gravitational force depends directly on Stot, 
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which is the fluctuation in the total matter density, while 
A describes the difference between the baryon fluctuations 
and 5tot ■ Before recombination, the baryons are dynamically 
strongly coupled to the photons while the dark matter fluc- 
tuations continue to grow independently. At lower redshifts 
the dominant contribution to the baryon fluctuation growth 
is the gravitational attraction t o the dark matter gravita- 
tional wells (e.g., see Figure 1 in iNaoz fc Barkanal (|2005T l). 

2.2 Large-Scale Structure: Small-fc Limit 

In the small-fc limit, the pressure terms can be neglected 
and equatio ns (|3| can be approxim ated as (following the 
derivation in lBarkana fc Loe b 200l) 

<W+2i3-<5 to t = l-ffo^r&ot , (4) 

A + 2HA = . 

Each of these equations has two independent solutions. As- 
suming that the universe is accurately described as EdS 
before reionization, these solutions are simple. The solu- 
tions for the Stot equation are growing and decaying modes 
-Dtot,i oc a and _Dtot,2 oc a~ 3 ^ 2 , while for the A equation 
we obtain Da.i oc 1 and Da, 2 oc a -1 / 2 . Moreover, with 
standard inflationary initial conditions, equations Q sat- 
isfy A(fc) oc Stot(k) at a given redshift. In other words, the 
relation between A and Stot is independent of wavenumber, 
for the range of redshifts and fc values that we are consider- 
ing here. We therefore find it useful to define 

r LSS = •5— , (5) 

Otot 

in terms of the solutions of equations ((4]) in the large-scale 
structure regime. The ratio tlss (which is negative) is inde- 
pendent of fc in this regime, and its magnitude decreases in 
time approximately oc 1/a, since A is roughly constant and 
5tot is dominated by the growing mode oc a. Figure [1] (top 
panel) shows |rLSs| as a function of redshift in the regime 
of large-scale structure. Although |rLSs| decreases as 2 — > 0, 
the initial difference between the dark matter and the baryon 
fluctuations has a large effect on the filtering mass even at 
2 < 10, as we show below. The small-fc regime (fc <C Uf) 
can be seen in Figure [T] (bottom panel), in terms of the 
filtering wavenumber fci? which is defined precisely in the 
following subsection. Note that remnants of the acoustic os- 
cillations in the photon-baryon plasma can be seen on the 
largest scales (k/kp ~ 10~ 4 , which corresponds to comoving 
fc ~ 0.01-0.1 Mpc -1 ), but these are much larger scales than 
those of high-redshift halos and the oscillations do not affect 
our results. 

A good numerical fit for n,ss in the redshift range of 
z — 7 ~ 150 is given by 

r-LSS = — H Vto + "3 (ft m ) , (6) 

where the form is motivated by the definition of tlss (equa- 
tion |[5j) and the EdS solutions of equations Q. The fitted 
coefficients are 

c*i(n m ) = 10" 4 x (-1.99 tf n + 2.41 n m + 0.21) , (7) 
a 2 (tl m ) = 10" 3 x (6.37 n 2 m -6.99 Q. m - 1.76) , 
a 3 (Qm) = 10" 2 x (-1.83 fi™ +2.4 ft m -0.54) , 




1+z 




k/k F 

Figure 1. The top panel shows |rxssl (equation J5]l) versus red- 
shift. The bottom panel shows the baryon to total fluctuation 
ratio versus k/kp for various redshifts. We consider, from bot- 
tom to top, z = 100, 65, 20, and 7. 

where the maximum residual error of the fit is 0.2%, over 
the range Q, m =0.25-0.4 in ACDM. 

2.3 The Filtering Scale 

In a ACDM universe, virialized CDM halos form on ex- 
tremely small scales at extremely early times. The minimum 
mass scale on which galaxies form within these halos is de- 
termined through a combination of two physical properties 
of the infalling gas. These are cooling and pressure, and each 
produces a characteristic minimum scale, so that the larger 
of the two scales becomes the dominant factor that deter- 
mines the minimum mass of star-forming halos. 

During the gravitational collapse, the gas potential en- 
ergy is transformed into kinetic and thermal energy through 
virialization and adiabatic compression. Unless the gas is 
able to dissipate its thermal energy with an effective cool- 
ing mechanism, the temperature will continue to rise, al- 
lowing the increasing pressure to halt the collapse. At very 
high redshifts, before metals were produced in supernovae 
and efficiently distributed in the intergalactic medium, the 
only cooling mechanism at temperatures below 10,000 K was 
cooling by molecular hydrogen, which itself is efficient only 
above a few hundred K. Thus, the first objects are expected 
to be fairly massive. Numerical calculations and simulations 
have shown that the minimum c ooling mass is ~ 10 5 M B at 
z ~ 100 and increases with time l|Tegmarklll997l ; lAbel et al.l 
120021 ; lFuller]|2000l ; lYoshidall2005l ). An object more massive 
than this minimum mass, after it goes through the virializa- 
tion process, will cool and collapse further, producing high- 
den sity clumps in which stars can form. 

lHaiman et al.l l|l997h showed that the UV flux needed 
to dissociate H2 in a collapsing environment is lower by more 
than two orders of magnitude than the flux that is necessary 
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to reionize the universe. Therefore, once stars reach a cer- 
tain abundance, the formation of stars through H2 cooling 
is suppressed, and atomic cooling becomes the only avail- 
able mechanism. This mechanism requires a much higher 
virial temperature (T vir > 10 4 K) which is associated with 
a halo mass of ~ 1O 8 M0. Therefore, in order to study the 
first generations of objects we must co nsider halo mass es of 
~ 10 5 - 1O 8 M (see also the review bv lBarkanal <|2006h ). 

On large scales (small wavenumbers) gravity dominates 
halo formation and pressure can be neglected. On small 
scales, on the other hand, the pressure dominates, and the 
baryon density fluctuations are suppressed compared to the 
dark matter fluctuations. The relative force balance at a 
given time can be characterized by the Jeans scale, which 
is the minimum scale on which a small perturbation will 
grow due to gravity overcoming the pressure gradient. If the 
gas has a uniform sound speed c 3 , then the comoving Jeans 
wavenumber is 



(8) 



Once the universe is matter-dominated, the Jeans scale is 
constant in time as long as T ~ 1/a, i.e., as long as the 
Compton scattering of the CMB with the residual free elec- 
trons after cosmic recombination keeps the gas temperature 
coupled to that of the CMB. At redshift z ~ 200 the gas 
temperature decouples from the CMB temperature and the 
Jeans scale decreases with time as the gas cools adiabati- 
cally. Any halo more massive than the Jeans mass can begin 
to collapse despite the pressure gradients. Figure [3] is a re- 
minder that the Jeans mass (dotted curve) is in the range 
10 4 - 1O 5 M0 during the formation of the earliest generations 
of galaxies. 

The Jeans mass is related only to the evolution of per- 
turbations at a given time. When the Jeans mass itself varies 
with time, the overall suppression of the growth of pertur- 
bations depends on a time-averaged Jeans mass. Following 
iGnedin fc Hull \ 1998), we define a "filtering" scale and use it 
to identify the largest scale on which the baryon fluctuations 
are substantial ly suppressed c o mpare d to those of the dark 
matter. While IGnedin fc Huil (| 19981 ) assumed a spatially- 
uniform sound speed and made a number of other approxi- 
mations (see Section f2. 5 [) . we calculate here the fitering scale 
obtained from the exact numerical solution of equations Q. 
In particular, our initial conditions at high redshift account 
for the coupling of the baryons and photons, i.e., 8b <C 5dm 
up to cosmic recombination. 

We define the filtering wavenumber Icf based on the 
scale at which the baryon-to-total fluctuation ratio drops 
substantially below its value on large scales. Thus, we ex- 
pand this ratio up to linear order in fc 2 , and write the ex- 
pansion in the following form: 



5 b 1 k 2 
x — = 1 ~ 7T + rLSS 

Otot K F 



(9) 



where tlss was defined in eq uation Our defini tion of fc^ 
is a generalization of that bv lGnedin fc Hull (|l99Sl ). who did 
not include the tlss term. To find fcf in our more general 
case we first write it in the form: 



**■(*) 



U(t) 



(10) 



where u(t) is to be determined. Substituting the expansion 
into equations ((!} and Q, we obtain an equation for u: 



u + 2Hu = / dl) 



1 k B f 



(<5tot + »"LSS<5tot + St) , 



(11) 



where we have neglected terms of higher order in fc 2 . We can 
solve this equation to find the parameter u: 



u(t) 



dt" 



dt'f d 



k B T(t') 



'Stot(t') + r hS s(t')5 tot (t') + ST(t')) 



(12) 



We have started the integral from the time of cosmic re- 
combination (t ICC ), since before recombination the contribu- 
tion of the baryon density and temperature fluctuations is 
negligible, i.e., the integrand essentially vanishes (Note that 
<*>tot + ^Lss^tot + St = 5b + St). As we discuss in Section [231 
this integral at high redshift (before cosmic reionization) 
gives a sig nificantly diffe r ent re sult from the approximate 
formula of IGnedin fc Huil (|l998h . 

Figure [T] (bottom panel) shows the baryon-to-total ra- 
tio versus fc/fci? at several different redshifts. The different 
values in the large-scale structure regime for different red- 
shifts are due to the different values of tlss (see top panel, 
Figure [TJ. The ratio drops on small scales. We have found a 
functional form that can be used to produce a good fit for 
the drop with wavenumber: 



= (1 + rxss 



+ TLSS 



(13) 



where x = fc 2 /kp, and n must be adjusted at each redshift. 
Defining 



V = 



1 + rLSS k F 
we can write the same fit as 



1 5 h 

1 + »"LSS <5tot 



y 



(14) 



(15) 



To second order, this gives: y ~ 1 — 77 + (n + l)r7 2 /(2n). 

We show in Figure [2] both the first and second-order 
approximations, as well as the full formula of equation (|15p . 
compared to the exact numerical results. We choose two red- 
shifts that bracket the interesting range, z — 100 (bottom 
panel) for the dark ages and z — 7 (top panel) for the latest 
possible beginning of reionization. From equation (|15[) we 
can see that the first-order approximation is independent of 
redshift. The figure also shows the exponenti al approxima- 
tion, g eneralized from the suggestion made bv lGnedin et al.l 
l|2003h for the post-reionization regime. In the universe prior 
to reionization, the exponential approximation y = exp(— rj) 
(which is also independent of redshift) becomes highly inac- 
curate at low redshifts. The figure makes clear that the drop 
of the fluctuations with k (or with 77) has a functional form 
that varies with redshift. With equation (|15|l we find that 
the power index at z = 100 is n = 23 (with behavior simi- 
lar to an exponential which would correspond to n — > 00), 
while at z = 7 it is n = 0.5. The values of n were found by 
matching the numerical result up to second order in r\. 

We define the filtering mass in terms of the filtering 
wavenumber using the traditional convention used to define 
the Jeans mass: 
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Figure 2. y versus r\ for two redshifts, z = 100 (top panel) and 
z = 7 (bottom panel). We compare the exact numerical results 
(dotted curves) to several approximations: the formula of equa- 
tion (1 1 5 I t (solid curves), the redshift-independent first-order ap- 
proximation (dot-dashed curves), the second-order approximation 
(short-dashed curves), and the redshift-independe nt exponential 
appro ximation (long-dashed curves) suggested by iGnedin et al.l 
1120031 ) for the post-reionization regime. 
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Figure 3. The filtering mass versus redshift. The exact calcula- 
tion (solid curve) , which can be fitted by equations II17I I and I I18H , 
is compared to a couple versions of the mean sound speed approx- 
imation (see Sections 12.51 and 12.61 1: where the integrals in equa- 
tion l|25| l are begun at time zero as written (long-dashed curve), 
or at cosmic recombination (short-dashed curve). Also shown is 
the Jeans mass (dotted curve), which is constant at z > 150. 



47T _ / 1 2n 



(16) 



Note that this relation, which we use consistently i n this 
paper, is one eighth of the definition in Gnedin (2000). The 
filtering mass essentially describes the largest mass scale on 
which pressure must be taken into account. Thus, in compar- 
ing different scenarios, those with higher gas temperatures 
will tend to have higher pressures, leading to higher values 
of the filtering mass. 

Figure [3] shows the evolution of the filtering mass with 
redshift. Since the baryon fluctuations are very small before 
cosmic recombination, the gas pressure (which depends on 
St) is small compared to gravity (which depends on dtot j 
see equations {T}). Thus, the filtering mass starts from low 
values and rises with time at z ~ 100. At lower redshifts the 
gas cools and the pressure drops. Therefore, even at z ~ 10 
the integral in equation (|12[1 receives a large contribution 
from much higher redshifts (z > 100). 

We have found a simple, accurate fit for the evolution of 
the filtering mass in the redshift range of z = 7 — 150. Using 
the notation Lm = \og(MF / Mq) and L z = log(l + z), the 
fit is of the form 

L M = /3i(n m )Ll+fa(n m )Ll+l33(n m )L z +l3 A {£l m ). (17) 



The fitted coefficients are 



/9i(fi m ) 

At(fi m ) 



-o.38 + o.4i n„ 



0.16 



(18) 



3.3 ntn - 3.38 Q. m + 1.15 , 
-9.64^ + 9.75 On -2.37 



= 9.8 ni 



10.68 fi m + 11.6 



where the maximum residual error of the fit is 0.2%, over 
the range ft m =0.25-0.4 in ACDM. 

2.4 Small Scale, Large-fc Limit 

In this limit the pressure dominates and makes 8b -C 
[<5tot , 5dm] , and therefore the first equation {1} becomes: 

<5 dm + 2H5dm — /dm^-ffo — J-<5dm • (19) 

Equation (|19|l has a simple analytical solution in the EdS 
regime: 5dm oc t a , where a = (— 1 ± yT+ 24/d m ) /6. In 
particular, the growing mode of the dark matter is reduced 
compared to the usual <5dm oc t 2 ' 3 solution, since while the 
baryons contribute to the cosmic expansion rate (in the 
Friedmann equation) they do not, on small scales, contribute 
to perturbation growth. Thus, since the dark matter fluctu- 
ations on large scales grow by a factor of ~ 300 between 
recombination and z = 10, then for /d m = 0.84 the to- 
tal growth of the dark matter fluctuations is reduced on 
small scales by a factor of ~ 1.8 relative to large scales. This 
linear calculation is limited to redshifts > 10, since even 
without reionization, by z = 10 the filtering scale becomes 
non- linear. 



2.5 Mean Sound Speed Approximation 

iNaoz fc Barkanal |2005l ) showed that the presence of spa- 
tial fluctuations in the sound speed modifies the calculation 
of perturbation growth significantly. Nevertheless, for com- 
pleteness and ease of comparison with previous results we 
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compare the above analysis to earlier, approximate calcula- 
tions. Thus, we proceed by applying a similar derivation as 
in the previous sections. In this approximation of a uniform 
sound speed, however, the evolution of the density fluctua- 
tions is described by a different set of coupled second order 
differential equations: 



5dm + 2H6 dm 
S b + 2H5 h 



Ho — 7T (/b<5b + /dm<$dm) 



(20) 



where c 2 = dp/dp is assumed to be spatially uniform (i.e., 
independent of k). With this assumption, the temperature 
fluctuations (as a function of k) are simply proportional at 
any given time to the gas density fluctuations: 



1 . 



&r_ __ c s 
5 h ksT/n 

This leads to different expressions for the filtering mass. 



(21) 



2.6 Mean Sound Speed Approximation: The 
Filtering Scale 

We again expand the ratio of the baryon fluctuations to 
dark matter fluctuatio ns in powers of k 2 . In this case we 
follow the derivation in lGnedin fc Huil l|l998ft and thus make 
several additional assumptions: that the baryon fraction is 
small (/b <C /dm), that the dark matter perturbation growth 
is dominated by the growing mode (<5dm oc D+), and that 
there is no initial difference between the baryon and the 
dark matter fluctuations. Although this derivation includes 
several different approximations, for simplicity we refer to 
it as the "mean sound speed approximation". With these 
assumptions, 



5dm 



= 1 



Writing this as 

5dm(t) 



(22) 



(23) 



and substituting into equations (|20p , we obtain an equation 
for the evolution of the parameter u Cs (analogous to equa- 
tion (JTTJ): 



ii c , + 2Hu c 



;D + (t) 



The solution (in analogy with equation (|12l) ~) is 



u Cs (t) 



dt" 
a 2 (t") 



dt'c 2 (t')D + (t') 



(24) 



(25) 



We note that the lower limit of the integral here is z — > oo, 
and not recombination as in equation (|12p . Before recom- 
bination, the coupling of the baryons with the radiation 
suppresses the baryon density fluctuations, but this is un- 
accounted for in this approximate calculation. Indeed, this 
formula implicitly assumes that the baryon perturbations 
grow like those of the dark matter, except for the effect of 
pressure. In particular, the filtering scale in this approxima- 
tion does not depend on the relative contributions of the 
baryons and the dark matter to the total matter density. 



In Figure |3j we show the filtering mass with this ap- 
proximation from iGnedin fc Huil fl998) . The correct filter- 
ing mass at z ~ 10 is substantially lower than predicted by 
the mean sound speed approximation. As noted before, this 
is mainly due to the fact that pressure forces depend on gra- 
dients in the gas density, and the baryon perturbations are 
reduced on all scales until they catch up with the dark mat- 
ter fluctuations (at z <C 100). In the approximate model, 
on the other hand, the initial difference between the baryon 
and the dark matter fluctuations is incorrectly neglected. 

In the Figure, we consider also the mean sound speed 
approximation where we start the integral in equation (|25[1 
at recombination (z — 1200). This should be more realistic 
than starting at time zero, since the baryon density fluctua- 
tion before recombination were negligible. This causes a rise 
in the filtering mass at z ~ 100, as in the correct calculation 
(i.e., the solid curve). However, the rise at high redshift is 
still much too fast, since this approximation still assumes 
that the baryons catch up with the dark matter fluctuations 
immediately after recombination. 

To summarize, the Figure shows that the difference be- 
tween the the correct calculation and the approximate one 
persists with time. This difference is remembered through 
the integrated effect of pressure, since at low redshift the 
Jeans mass is lower and thus the pressure is lower as well. 
Therefore, the high redshifts contribute most to the inte- 
grated pressure, and thus even though the difference be- 
tween the baryon and the dark matter fluctuations declines 
with time (e.g., see Figure [T] top panel), the system still 
remembers the initial difference. 



2.7 Gas Fraction 

One useful application of the filtering mass is to the esti- 
mation of the fraction of gas inside halos. iGnedinl (|2000| . 
his equation (8)) estimated the mean baryonic mass M g in 
halos of total mass M to t using a formula fitted to his post- 
reionization simulation: 

/bM tot 



Mg(M t0 t,t) 



(26) 



[1 + (2V3 _ l)M c (t)/M tot f 

where M c (t) is a characteristic halo mass that corresponds 
to a gas fraction of 50% of the cosmic baryon fraction. It 
is natural to expect a close relation between the character- 
istic halo mass and the filtering mass, since the gas frac- 
tion in a collapsing halo reflects the amount of gas that was 
able to accumulate in the central, collapsing region, during 
the entire extended collapse process. In particular, if the 
Jeans mass changes suddenly, this does not immediately af- 
fect then-collapsing halos. A change of pressure immediately 
begins to affect gas motions (through the pressure-gradient 
force), but has only a gradual, time-integrated effect on the 
overall amount of gas in a given region. 

As mentioned before, our definition of the filtering mass 
(equatio n (I16|> ) is one eighth of the previous definition of 
IGnedinl 112000ft. T hus the characteristic mass that matches 
the lGnedinl (|2000h simulations corresponds to about 8 times 
our filtering mass. Assuming that this is true also prior to 
reionization, we can evaluate the total gas fraction in halos 
as a function of redshift as 



F a {z) 



fs T (M tot ,z] 



M a {M totl t) 
M to ± 



dS 



(27) 
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where S = a 2 ( M,z) is the variance and fsr is the 
ISheth fc Tormenl l|200ll ) function for the fraction of mass 
associated with halos of mass M (explicitly given in equa- 
tion (|33[) . below). 

The total fraction of gas that is in halos is shown in Fig- 
ure^ for our correct calculation of the filtering mass as well 
as for the previous calculation (which we have referred to as 
the mean sound speed approximation). We also compare to 
the fraction of gas in halos above the minimum H2 cooling 
mass, and to the fraction above the minimum atomic cooling 
mass. In the correct calculation, a significant part (10-50%) 
of the total gas in halos arises from halos that are below 
the characteristic mass. The prediction of the gas fraction 
in halos in our correct calculation is higher than that based 
on the previous approximation, by a factor > 2 at high red- 
shifts and still by 10% at redshift 7. In this redshift range, 
around half the gas in halos is in potentially star-forming 
halos (i.e., those with efficient H2 cooling), and the rest is 
in gas minihalos. In particular, this means that the small- 
est star-forming halos (i.e., those with a mass equal to the 
minimum H2 cooling mass) are moderately affected by pres- 
sure, and have their gas content reduced to around half the 
cosmic baryon ratio. 

The importance of the pressure in halos of mass equal 
to the minimum H2 cooling mass is illustrated in Figure 
We consider the improved calculation (solid curve) and the 
mean sound speed approximation (short-dashed curve). This 
Figure shows that for the very first stars the previous cal- 
culation underestimates the gas fraction in these halos by 
more than an order of magnitude. E.g., at 1 + z — 66 (the 
formation of the first star) the improved calculation predicts 
that the gas is about 35% of the cosmic baryon fraction in 
halos with mass equal to the H2 cooling mass, while the 
mean sound speed approximation predicts only 1.3%. Thus, 
we conclude that the effect of pressure on the very first stars 
is only moderate, unlike the result suggested by the mean 
sound speed approximation. The discrepancy decreases with 
the redshift. However, at 1 + z = 20 the previous calculation 
still underestimates the gas fraction by a factor of 2, and 
even at 1 + z — 8 the previous calculation underestimates 
the gas fraction by about 10% compared with the improved 
calculation. 



3 FORMATION OF NON-LINEAR OBJECTS 

The small amplitude density fluctuations probed by the 
CMB grow over time as described by equations ([l]) and ([2]), 
as long as the perturbations are linear. However, when S 
in some region becomes of order unity, the full non-linear 
problem must be considered. The standard calculation that 
describes the formation of spherical non-linear objects was 
done for a dark matter halo only, as was explained in Sec- 
tion \T\ Here, we generalize the calculation to the high- 
redshift regime, including the gravitational effects of the 
baryons and the radiation. 

The initial conditions in N-body simulations of the first 
galaxies are set long after the recombination epoch (usually 
at z ~ 200 or later). Halos that collapse early (z ~ 20 and 
earlier) thus cannot arise from fluctuations that are linear 
at the beginning of the simulation, regardless of the size 
of the corresponding perturbed region. Now, perturbations 




1+z 



Figure 4. Top panel: Gas fraction in halos versus redshift. We 
show our correct calculation (solid curve) and the previous cal- 
culation using the filtering mass in the mean sound speed ap- 
proximation (dotted curve). We also compare to the fraction of 
gas in halos above the H2 cooling mass (long-dashed curve), and 
the fraction above the atomic cooling mass (short-dashed curve). 
Bottom panel: Ratio of gas fractions in different cases. We con- 
sider the gas fraction in halos above the characteristic mass (solid 
curve), the gas fraction in halos above the H2 cooling mass (long- 
dashed curve), and the gas fraction in the mean sound speed 
approximation (dotted curve); each is divided by the total gas 
fraction in halos in our correct calculation. 



that correspond to galaxies or clusters start from a comov- 
ing length below 100 Mpc and are thus expected to have 
entered the horizon at some time t ente r in the radiation 
dominated era. The perturbations then do not grow signifi- 
cantly until equality (z ~ 3000). If we follow the evolution 
of a collapsing halo, then the linear <5 grows approximately 
oc a after equality. Thus a region that collapses to a halo at 
z ~ 100 must have entered equality already in the weakly 
non-linear regime (5 ~ 1.686/30 ~ 5%); since growth was 
slow before then, even when the halo entered the horizon in 
the radiation dominated universe, the perturbation was not 
extremely small. In practical applications, we find that the 
largest 8 gets at hori zon crossing is ~ 10~ 3 , for halo s hosting 
the very first stars ( Naoz. Noter fc Barkanall2006l ). so that 
non-linear corrections are always small. For completeness, 
however, we first consider the behavior of spherical fluctua- 
tions outside the horizon. 



3.1 Fluctuation Growth Outside the Horizon 

We consider a spherical top-hat fluctuation bigger than the 
Hubble radius at an early cosmic time. Since the fluctuation 
is outside the horizon, we cannot use Newtonian perturba- 
tion theory, and must apply General Relativity. 

We consider a spherical overdensity with a uniform den- 
sity p = p (1 + S). Birkhoff's theorem implies that the Fried- 
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Figure 5. The gas fraction in halos of mass equal to the mini- 
mum H2 cooling mass compared to the cosmic fraction. We plot 
the improved calculation (solid curve) and the mean sound speed 
approximation (short-dashed curve). In the latter case we have 
used the previous estimate, starting the integrals in eq. (12 5 I t at 
t = 0. 



mann equation is the correct solution within the spherical 
region (which has a curvature k): 



H 



sph 



{1 + 5) + ^(i + 5 )i + 



a sph 



.(28) 



where H sp h and a sp h are the Hubble constant and expansion 
parameter, respectively, associated with the evolution of the 
perturbed region, while a describes the cosmic expansion of 
the mean universe. The 4/3 power is a result of the differ- 
ence between the non-relativistic (dark matter and baryon) 
density fluctuation 8 and the radiation density fluctuation. 

Just as in the regular top-hat collapse, we wish to com- 
pare the exact non- linear evolution (given by equation (1281) ) 
to the linearly-extrapolated evolution. Comparing the above 
non-linear equation to the evolution of the background uni- 
verse and linearizing, we obtain: 



2Ha8 = I 2 S 



k 



(29) 



where I 2 = H 2 [fl m /a 3 + (4/3)fi r /a 4 ], and a = 0.5 I 2 /H; at 
high redshift, in the radiation-dominated regime, a « —0.25 
and then 8 oc a 2 (corresp o nding to the synchron ous gauge; 
iMa fc Bertschingerl (|l995l ) ; iPadmanabhanl (|2002l )). We start 
early enough so that 8 — 1 x 10~ 4 initially, and assume the 
growing mode for the initial perturbation. Solving numer- 
ically the non-linear and linear equations (equations Q280 
and (J29J) ) , until the fluctuation enters the horizon, yields 
the initial value of the fluctuation entering the horizon in 
the radiation dominated universe. As noted above, in prac- 
tice, linear initial conditions at horizon crossing would suffice 
for the parameter space we consider in this paper. 



3.2 Fluctuation Growth Inside the Horizon 

In the Newtonian regime the non-linear growth is described 
by the Newtonian equation (or, more precisely, the equa- 
tion for the acceleration that results from the Einstein equa- 
tions): 

' = (30) 

where the rest stands for all matter that does not participate 
in the collapse, and thus only contributes to the expansion 
of the universe. We define rd m and Vb to be physical radii 
that enclose a fixed mass of dark matter and of baryons, 
respectively, assuming a tophat perturbation in each. Then 
we obtain two coupled non-linear equations of motion: 



ran 



-1 4ttG 



r 2 3 

' dm 



fdm (Pdm + pb) + -HofiArdn 



8ttG 



-1 4ttG 



Tb (Pdm + Pb) + H Q A rb 



8ttG 



- pr fdm 



p r r b ■ (31) 



We have assumed that the radiation is kept smooth by 
its own pressure and does not participate in the collapse; 
the factor 87r/3 is the result of inserting P r = p r /3 in 
equation (|30[) . Since the time when the fluctuation en- 
ters the horizon is substantially early in the radiation- 
dominated universe, the baryon-photon coupling yields 
Sb, S-, <C <5dm,5tot initially. We calculate separately the lin- 
ear and non-linear growth of the fluctuations. The resulting 
critical (linear) overdensity at the time of collapse is shown 
in Figure [6] For the cosmological parameters that we use in 
this paper, we find that S c is essentially independent of M, 
and is lower than the EdS value by ~ 1% x (1 + z)/20 in 
the range of z — 9 — 100. When dealing with very rare ha- 
los, even a change of a few percent in S c can change the halo 
abundance at a given redshift by over an order of magnitude 
(see Figure 2 of lNaoz. Noter fc Barkanal l|2006l )). The results 
are insensitive to the cosmological parameters: also shown 
in Figure HJ] are the resu lts for the s et of other parameters 
(Section 1) specified by IViel et ail (|2006l ). which differ by 
1-a from our standard WMAP parameters; the results are 
almost identical in the EdS regime (1 + 2 = 2 — 10), with 
only a < 0.2% difference at high redshift (and even less at 
z < 1). 

The effect of various physical ingredients on 8 C can be 
illustrated using the different cases shown in the upper panel 
of Figure [S] In general, collapse is most efficient in the EdS 
case where all the matter participates in the collapse (re- 
sulting in 8 C — 1.686); any smooth component that does not 
collapse (a cosmological constant at low redshift, radiation 
and baryons at high redshift) reduces the collapse efficiency 
since only the dark matter component takes part in the col- 
lapse throughout. Now, any reduction in the matter fraction 
that collapses depresses the linear evolution of the density 
perturbation more strongly, while the non-linear perturba- 
tion is larger and is thus less affected by the components 
that do not help in the collapse. Therefore the linear per- 
turbation reaches a lower value of 8 C when the non-linear 
perturbation collapses. As shown in the upper panel of Fig- 
ure |6j the extended period at high redshift when the baryon 
perturbations remain suppressed is the main cause of the 
reduction of the value of <5 C , but the contribution of the pho- 
tons to the expansion of the universe also makes a significant 
contribution. 
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Figure 6. Critical overdensity versus collapse redshift. We com- 
pare 8 C for the full calculation (solid curves) to the results with 
other assumptions. The top panel compares to several cases that 
include only some of the physical ingredients that affect the spher- 
ical collapse calculation. We show the results of previous cal- 
culations that did not properly include the baryons and pho- 
tons at high redshift (dotted curve), the results if the baryons 
are treated correctly but the contribution of the uniform radi- 
ation background is neglected (long-dashed curve), and the re- 
sults if the radiation is included but the baryon-photon cou- 
pling is neglected (i.e., the baryons are treated like dark matter: 
short-dashed cur ve). The bottom panel shows results also for the 
IViel et all ([2006) set of parameters (dashed curve). 



3.3 The Mass Function and Halo Abundance 

In addition to characterizing the properties of individual ha- 
los, it is important for any model of structure formation to 
predict also the abundance of halos. We calculate in this 
section the number density of halos as a function of mass, 
at any redshift. We start with the simple analytical model of 
iPress fc Schechterl l| 19741) , which is based on Gaussian ran- 
dom fields, linear growth and spherical collapse. It predicts 
that the abundance of halos depends on mass and redshift 
through the two functions a(M, z) and S c (z), where a 2 (M, z) 
is the variance (calculated from the power spectrum) as a 
function of halo mass at z, and S c (z) is the critical collapse 
overdensity from Section \3. 2 1 

In Figure [7] we show the number of a that a fluctuation 
must be in order for it to collapse, at various redshifts from 
the first star until the present. For example, a 1O 5 M0 halo 
that collapses and forms the first star is a 9.5-ct fluctuation 
according to this criterion (but see below). Halos of mass 
10 5 Mq become l-cr fluctuations only at z ~ 5, while today 
the l-cr collapsing fluctuations are 10 13 Mq group halos. 

The growth of the density fluctuations in our correct 
calculation can be described by an effective growth factor 
which depends on the mass scale: D e g = a(M,z)/a(M,z = 
0). We compare this e ffective growth factor to the traditional 
growing mode D(z) jPeeblesl Il980l ) in Figure [8] where we 



Figure 7. Rarity of fluctuations that produce halos at various 
redshifts. We show the size of a fluctuation that produces a col- 
lapsed halo on the scale M at redshift z, in terms of the typi- 
cal fluctuation level (i.e., measured as a number of cr). This in- 
cludes <j(M, z) from the cor rect calculation of the power spec- 
trum IINaoz &; B arkana 2005), and S c (z) from the correct spher- 
ical collapse calculation (Section l3.2l l, We consider, from bottom 
to top, 2 = 0, 1.2, 6.5, 11, 47 and 66. Redshift 66 corresponds 
to the formation of the first star, while 2 = 47 corresponds to 
the redshift of the second generation of stars, i.e., the first col- 
lapse via atomic cooling. Redshift 11 corresponds to the first 
halo as massive as that of the Milky Way and 2 = 1.2 is associ- 
ated ^with_Uiefc£mationoftliefirst cluster as massive as Coma 
fsee lNaoz, Noter fc Barkanal j2006h l. We also show the results at 
2 = 6.5 associated with observations of the most distant quasars 
(and perhaps with the end of reionization). 



show D c ff /D—l as a function of mass, illustrated at the same 
redshifts as in Figure[7] The source of the difference is in the 
smoothness of the radiation and baryons at high redshifts; in 
addition, pressure suppresses growth on small mass scales. 
These effects make D e g larger at high redshift compared to 
its value today (since there has been less growth). 

The comoving number density of halos of mass M at 
redshift z is 



dn _ po „ 
dM ~ M /ST 



<1S 



dM 



(32) 



were we have used the ISheth fc Tormenl l|200ll ) mass func- 
tion that fits simulations, and includes non-spherical effects 
on the collapse. The function fsr is the fraction of mass 
associated with halos of mass M: 



1 + 



(a'u 2 )"' 



exp 



(33) 



where v = 8 c /\fS. We use best-f it parameters a — 0.75 and 
q' = 0.3 (|Sheth fc Tormenll2002i ), and ensure normalization 
to unity by taking A' — 0.322. We apply this formula with 
Sc(z) and a 2 (M,z) as the arguments. The lSheth fc Tormenl 
|200il ) mass function makes it easier for rare fluctuations 
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Figure 8. The fractional difference between the effective growth 
factor of our correct calculation, cr(M, z)/cr(M, z = 0), and the 
traditional growth factor D(z), as a function of mass at various 
redshifts. We consider, from bottom to top, z = 0, f.2, 6.5, ff , 47 
and 66, as in Figure [71 

to collapse compared to lPress fe Schechterl (Il974l ). This, for 
instance, makes the first star-forming halo effectively only 
an 8.3-cr density fluctuation on the mass scale of 10 s M , m 
terms of the total cosmic mass fraction contained in halos 
above this mass. 

The cumulative comoving number density of halos n(> 
M min ) is given by 

n(>M^)=J M ^ d M. (34) 

In Figure [9]we plot the cumulative number density of halos. 
In addition to our standard WMAP cosmological para me- 
ter set, we have also compared to the lViel et alJ (|2006h set 
of parameters (Section 1). We obtained for these parame- 
ters that the first star formed at z = 60.0 (instead of 65.8) 
and the first Coma size cluster formed at z — 1.06 (instead 
of 1.24). The difference between the two sets of cosmologi- 
cal parameter arises mainly from the difference in as. The 
difference in the values of S c is negligible (see Figure [6]) . 



4 CONCLUSIONS 

We have calculated the filtering mass correctly at high red- 
shift and compared it to previous estimates. We have found 
that at high redshift the filtering mass is lower by about 
an order of magnitude compared to previous calculations. 
The difference declines with time but remains a factor of 
~ 3 even at z = 7 (Figure [3}. Our calculation predicts a 
lower filtering mass because it includes the initial difference 
between the dark matter and baryon fluctuations, due to 
the baryon-photon coupling, and this lowers the pressure. 
This means that in contrast with the previous prediction, 



Figure 9. The halo abundance as a function of the halo mass at 
different redshifts. From top to bottom (at the high-mass end): 
z = 0,1.2,6.5,11,47 and finally z = 66 which corresponds t o 
the formation of the first star llNaoz. Noter & Barkanal l2006h . 
We consider our standard W MAP cosm o logica l parameters (solid 
curves) , and compare to the IViel et all l l2006t) set of parameters 
(dotted curves). 



pressure has only a moderate effect on the formation of the 
earliest luminous objects. 

Before recombination the baryon fluctuations were sup- 
pressed compared to the dark matter fluctuations due to 
tight coupling with the radiation. The filtering mass rises 
with time at high redshift because after recombination the 
baryon pressure gradients start to increase. However, after 
the baryon temperature decouples from the CMB tempera- 
ture the gas cools adiabatically, the Jeans mass drops, and 
eventually (at lower redshifts) the filtering mass drops as 
well (see Figure [3} . The high pressure at very high redshift 
still contributes significantly to the filtering mass at redshift 
below 10. The delay in the drop of the filtering mass com- 
pared to the Jeans mass is a signature of the continuing 
contribution of the memory of the pressure from very high 
redshifts. We have found numerical fits to the difference be- 
tween the baryon and the total density fluctuations on large 
scales (equations (0); we have also fit the filtering mass as 
a function of fi m and z (equations (|17|) and (|18[0 , 

Using the new prediction of the filtering mass we have 
shown that at high redshift there is more gas in halos than 
in previous estimates. This difference is as large as a factor 
of 2, but remains ~ 10% even at redshift 7. Around half of 
the gas is in halos with efficient H2 cooling, and the rest is 
in gas minihalos. The previous calculation also suggests a 
greatly reduced gas fraction in the halo that hosts the first 
star, while we find only a moderate effect (Figure [5}. 

In addition, we have computed the evolution of linear 
as well as non-linear spherical overdensities, outside and in- 
side the horizon, in a ACDM universe. Previous analyses 
showed that the cosmological constant contribution to the 
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expansion of the universe results in a drop of the value of 
S c by ~ 0.6% today (e.g.. iLahav et al] [l991). This makes a 
significant difference in the abundance of clusters. Consid- 
ering structure formation at high redshift, we investigated 
the effect on 8 C and on the halo abundance of the contri- 
bution of radiation to the expansion of the universe, and of 
the contribution of the baryons to the collapsing halo, given 
their different initial conditions compared to the dark mat- 
ter. We have found that there is a 3% change in the value of 
the overdensity S c at z ~ 60 (Figure [6]). This changed value 
translates to a much larger difference in the halo abundance 
at high redshift (Figure [5]); these differences decline at low 
redshift. The large difference in the halo abundance is the 
result of the difference in the exponent in the mass function 
(equation (|33p L and shows that even small effects on halo 
formation can in some cases be very important. 
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